Gain-induced trapping of microcavity exciton polariton condensates 
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We have performed real and momentum space spectroscopy of exciton polariton condensates in 
a GaAs-based microcavity under non-resonant excitation with an intensity stabilized laser. An 
effective trapping mechanism is revealed, which is due to the stimulated scattering gain inside the 
finite excitation spot combined with the short lifetime. We observe several quantized modes while 
the lowest state shows Heisenberg-limited real and momentum space distributions. The experimental 
findings are qualitatively reproduced by an open dissipative Gross-Pitaevskii equation model. 
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It is well known that an infinite two-dimensional (2D) 
system cannot establish true long-range order above zero 
temperature because of phase fluctuations But 
when a trapping mechanism is employed, so that the sys- 
tem size is finite, a crossover temperature exists below 
which there are sizeable correlations across the trap. For 
example, a 2D system of ideal bosons in a trap is ex- 
pected to show Bose Einstein condensation (BEC) below 
a temperature that depends on the trap size and shape 

i- 

Microcavity exciton polaritons Q is a 2D system that 
behaves as bosonic at low densities. Although the parti- 
cle lifetime is short, limited by the photon lifetime inside 
the cavity, so that true thermal equilibrium cannot be 
established, a dynamical condensation phenomenon sim- 
ilar to BEC is observed [HI, ^ . Free electron- hole pairs 
are optically injected, which then relax towards the lower 
polariton (LP) branch. Above a threshold laser pump- 
ing power the scattering rate towards one or a few final 
states is significantly enhanced due to the bosonic final 
state stimulation, and condensation occurs. It has been 
proposed that external trapping such as photonic disor- 
der 0, stress trap [H , micropillar patterning Q , or metal 
film [10| is essential to stabilize a spatial mode. 

Here, we show that there is an implicit confine- 
ment mechanism due to the finite lifetime and the spa- 
cially confined gain inside the finite excitation spot size. 
Namely, reservoir polaritons cannot travel far from the 
point where they were created and final state stimulation 
only occurs inside the laser excitation spot area where 
polariton density is sufficient. Under suitable excitation 
conditions, we can observe the quantized states of this ef- 
fective trap in real as well as in momentum space. Only 
the lowest state is perturbed by the disorder potential 
and it is further confined by the latter. We have observed 
Heisenberg-limited real and momentum space distribu- 
tions for this state, which confirms that phase fluctua- 
tions are negligible. We use a simplified one-dimensional 
model, employing a generalized Gross-Pitaevskii equa- 



tion with loss and gain terms, to qualitatively reproduce 
our experimental results. 

The sample is the same as in our recent experiments 
111 Il2j , and it consists of an AlAs | cavity sandwiched 



between two distributed Bragg reflector (DBR) mirrors. 
The upper and lower mirrors are made of 16 and 20 
pairs respectively of AlAs and Gao.8Alo.2As. 3 stacks 
of 4 GaAs quantum wells (QW's) are grown at the cen- 
tral three antinodes of the cavity. The measurements 
reported here are taken from a spot on the sample with 
photon-exciton detuning S = —2meV^ while the Rabi 
splitting is 2hQ.Rahi = ^^meV. The sample is kept at a 
temperature of 7 — on the cold finger of a He flow 
cryostat. The system is pumped with a single- frequency 
Ti- Sapphire ring laser focused on a flat spot of diame- 
ter ~ 25/im (as measured for low pumping power) from 
the direction normal to the sample [li^. The laser is lin- 
early polarized, and the orthogonal linear polarization is 
detected. Luminescence above threshold is almost unpo- 
larized, namely the degree of linear polarization is < 0.15. 
This is attributed to the small ground state linear polar- 
ization splitting, typically ~ SO/ieF in our sample [ll| . 
We expect that the condensate randomly chooses one di- 
rection of linear polarization in every realization, there- 
fore luminescence appears almost unpolarized after time- 
averaged detection. The laser wavelength is tuned to the 
first reflectivity minimum of the cavity above the stop 
band and it is modulated into pulses with O.Sms dura- 
tion and lOOi^z repetition rate to prevent sample heat- 
ing. Luminescence is collected through an objective lens 
with numerical aperture NA = 0.55. The spectroscopy 
setup is the same as in [ll| and allows us to perform 
near field (NF - real space) and far field (FF - momen- 
tum space) imaging and spectroscopy. That is, we can 
measure energy-resolved luminescence as a function of 
position or of in-plane momentum. Our resolution is Ifim 
for NF spectroscopy and 0.05/im~^ for FF spectroscopy. 

Above a threshold pumping power of 70m W^, the phase 
transition from normal to condensation states is observed 
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FIG. 1: (color online) Evidence for polariton condensation, 
(a) Power emitted between ±4° as a function of pumping 
power. A non-linear increase is evident at a threshold power 
of 70mW. (b) Normalized spectra as a function of power. 
Above threshold, several narrow peaks appear. Inset: Energy 
shift of the main spectral peak as a function of power. The 
linear blue shift below threshold is followed by a logarithmic 
shift above threshold. 
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FIG. 2: (color online) (a) Real space and (b) momentum space 
spectra measured above threshold (at 115mW). 



(Fig[T]): the power emitted at small angles increases non- 
linearly, spectral peaks with reduced linewidth appear 
and the energy blue shift changes from a linear depen- 
dence on the pumping power to a logarithmic dependence 
[9[. Fig. [2] shows real space and momentum space spec- 
tra measured above threshold. Six discrete modes ap- 
pear whose linewidth is resolution-limited. Our spectral 
resolution has HWHM ~ 25/iey, corresponding to co- 
herence time h/ HWHM = 26ps. 

In order to move the real or momentum space image 
along one direction, we move the last lens of our optical 
system (the one in front of the spectrometer) along the 
same direction perpendicular to the optical axis. The 
translation is also perpendicular to the direction of the 
spectrometer slit. This way, we can record spectra along 
parallel lines in real or momentum space resulting in a 
2D grid of spectrally-resolved luminescence. We can then 
integrate around the energy of each mode and generate 
2D images of the different modes in real (Fig. [3fa)) as 
well as in momentum (Fig. HJ^a-f)) space. 

Inside a circular trap with infinite barrier height, the 
solutions of the Schrodinger equation are expressed in 
terms of Bessel functions of the first kind 

V'„^„(r,,^) ae^™'^J™(fc„r), fc„ = (1) 

where Rq is the trap radius and jmn are the roots of 
the Bessel function Jm(^) for ^ = 1,2,3 — The FF 
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FIG. 3: (color online) (a) Experimental 2D real space image 
of the lowest trapped state with quantum numbers (m,n) = 
(0,1). (b-g) ID mode profiles, blue: experiment, red: trapped 
states of an infinite circular quantum well of radius 8fim. 
In (b) we have integrated along the ^/-axis, and also fitted 
the experimental profile with a gaussian distribution (black 
dashed line), (c-g) are simple cross sections along the x-axis 
ai y — 0/jm. The quantum numbers {m,n) are (b) (0, 1), (c) 
(1,1), (d) (2,1), (e) (0,2), (f) (1,2), and (g) (0,3). 
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FIG. 4: (color online) (a-f) Experimental 2D momentum 
space images of trapped states. The quantum numbers (m, n) 
are (a) (0,1), (b) (1,1), (c) (2,1), (d) (0,2), (e) (1,2), and 
(f) (0,3). (g-1) ID mode profiles, blue: experiment, red: 
trapped states of an infinite circular quantum well of radius 
8/im. In (g) we have integrated along the /cy-axis, and also 
fitted the experimental profile with a gaussian distribution 
(black dashed line), (h-1) are simple cross sections along 
ky = 0/im~"^. The quantum numbers {m,n) are (g) (0,1), 
(h) (1,1), (i) (2,1), (j) (0,2), (k) (1,2), and (1) (0,3). 



wavefunctions follow from the Fourier transform of the 
above expression. Based on this, we can assign quantum 
numbers (m, n) to the six observed modes as shown in 

Figs. [312 

The real space images show the characteristics of 
the above wavefunctions. The comparison is shown in 
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Fig. [3jb-g) for a trap of radius 8/im. This radius is 
shorter than the spot radius measured below threshold 
(~ 12.5/im), as the spot shrinks above threshold. This 
is because the lifetime of the excitations we create drops 
from ~ 200^5 below threshold to a few ps above thresh- 
old, as measured for our sample and excitation condi- 
tions, so that they do not have enough time to diffuse. 
This is due to the bosonic final state stimulation effect, 
namely the relaxation rate is proportional to the pop- 
ulation of the final state, so it is accelerated above the 
condensation threshold [13|- For example, a polar iton 
of mass 6 x 10~^me and momentum kx = 0.5/im~^ can 
only move by 2/im in 2ps. The inability of polaritons to 
escape from the pumping spot creates a finite-size gain 
region that traps condensates inside it. This is a sim- 
ilar mechanism to gain guiding in semiconductor laser 
structures 

The same model works for the momentum space im- 
ages, as shown in FigHJ^g-l). We note that the harmonic 
potential model could not explain the energy splitting be- 
tween modes (2,1) and (0,2) (3rd and 4th modes) which 
is clear in Fig. [2jb). Also, the modes (m,n) = (3,1), 
(4,1), (2,2), which should appear among the observed 
modes, are missing since their large angular momentum 
makes polariton scattering into these states inefficient. 

The sample disorder potential follows a gaussian distri- 
bution with standard deviation < O.lmeV [12| which is 
considerably smaller than the kinetic energy of all modes 
except the lowest one, as shown on the measured spectra 
of Fig. [21 Indeed, only the lowest mode is perturbed 
by the disorder potential and it is further confined. Fig. 
[31(a) shows that the lowest mode has migrated away from 
the center of the beam, possibly towards a local minimum 
of the disorder potential. In Fig. ^h) we plot the mea- 
sured profile of this mode along the X axis . where we 
have integrated Fig [31(a) along the axis and shifted 
the result to center it around x = 0/im (Instead, Figs. 
[3l(c-g) are simple cross sections along y = Ofim). The 
observed distribution is narrower than the expected one 
based on the infinite quantum well model. We can fit 
it with a gaussian (dashed line) with standard deviation 
(Jx = 2.16/im. The same effect is evident in the mo- 
mentum space distribution of the lowest mode (Fig [41(g)) 
which is wider than the model prediction. Fitting with a 
gaussian gives a standard deviation of = 0. 283 fim~^ . 
Taking the product cix x (Jk^ = 0.61, which is very close 
to the Heisenberg limit of 0.5. 

The stabilized intensity of the excitation laser is criti- 



Ti-Sapphire laser operated in the CW mode. As the laser 
cavity is designed for generating psec pulses, it supports 
a large number of modes with frequencies inside a 
wide window, so we expect power fluctuations at the ps 
timescale. Power fluctuations cause fluctuations in the 
polariton density and the interaction energy is modulated 
accordingly. This is the main decoherence mechanism 



[15|. Indeed, the coherence time of the polariton con- 
densate measured with a Michelson interferometer was 
Tco/i < 5ps in that case. In the present experiment, the 
laser cavity is designed for single-mode operation and a 
double etalon is placed inside it to suppress all but one 
mode. The coherence time of the lowest observed polari- 
ton mode is Tcoh > 26ps as evidenced from the resolution- 
limitted linewidth in Fig [21(b). Measurement with the 
Michelson interferometer gives the result Tcoh ^ 20ps, 
but the real coherence time is expected to be longer, as in 
this measurement we could not fully suppress the higher 
modes. This is consistent with the results reported in 
[l6j on a different sample. 

We see that pumping with such a laser source cre- 
ates the appropriate conditions for a steady state con- 
densate with long coherence time. This is why only the 
eigenmodes of the gain-induced trap are visible above 
threshold. On the other hand, when the laser has power 
fluctuations on short timescale, it introduces an energy 
uncertainty as explained above. In this case, polaritons 
condense into a mixed state composed of different modes 
closely spaced in energy. In particular, phase fluctuations 
are prominent [HI . The situation is more complicated in 
highly disordered samples [l^ , in which case the disorder 
potential plays a big role in determining in which states 
polaritons condense. 

The non-thermal population distribution among the 
different modes is due to the non-equilibrium character of 
polariton condensation and is determined by the pump- 
ing and loss rate of every state. In particular, when we 
change the beam profile from a Il-shape to an M-shape, 
so that the pumping density is higher far from the cen- 
ter, the lowest modes are suppressed. The Heisenberg- 
limited NF and FF distributions of the lowest mode con- 
firm that this state is highly coherent, without substan- 
tial phase fluctuations that would introduce energy un- 
certainty. The interference patterns observed in the FF 
distributions of the higher modes confirms the coherent 
character of these states, as well. 

We employ a one-dimensional theoretical model sim- 
ilar to the one used in 0. It cosists of two coupled 
equations, the open dissipative Gross-Pitaevskii equation 
for the condensate order parameter ip{x^t) and the rate 
equation for the reservoir population URix^t) 



ih 



dil){x^ t) 
dt 



+ Vext{x) - —[iC - RnR{x,t)] 



cal for our observations. In [12|, we used a mode-locked dnR{x,t) 



dt 



gc\i^{x,t)f ^ gRnR{x,t)}ilj{x,t), (2) 

Plas{x,t) - jRnR{x,t) 

RnR{x,t)\^{x,t)\\ (3) 



Vext{x) is the external disorder potential and is the 
condensate decay rate. For the stimulated scattering rate 
from the reservoir to the condensate RxnR{x^ t), we use a 
linear dependence on the reservoir population, gc and gR 
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FIG. 5: (color online) Typical results of our theoretical model 
(a) The disorter potential Vext{x) and (b) the laser pumping 
rate profile Pias{x) we considered, (c-d) Real space spectra for 
different pumping powers, 2 and 10 times respectively above 
threshold. Only the lowest mode is confined by the disorder 
potential. 



complete calculation in a two-dimensional grid, and is be- 
yond the scope of the paper. 

We have presented evidence for an inherent effective 
trapping mechanism of exciton polar itons, which is due to 
the finite size of the excitation spot along with the short 
lifetime. Under excitation with an intensity-stabilized 
laser, several quantized modes of the eflFective trap spon- 
taneously appear, whose kinetic energies are higher than 
the disorder potential fluctuations. We observed this 
quantized mode structure both in real and in momen- 
tum space spectra and showed that the real and momen- 
tum space distributions of the lowest mode are Heisen- 
berg limited. A one-dimensional model based on the 
open dissipative Gross-Pitaevskii equation with loss and 
gain qualitatively reproduces our experimental flndings. 
This trapping mechanism should be present to all sam- 
ples studied so far. 
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are the condensate-condensate and condensate-reservoir 
interaction constants. The reservoir density riRix^t) is 
controlled by the laser pumping rate Pias{x^t) and the 
reservoir decay rate 7^^. In the time-dependent simula- 
tion, ilj{x^t) is a superposition of terms oscillating at dif- 
ferent frequencies, which correspond to different modes. 

For our simulations, we set the interaction parame- 
ters gc and qr equal to zero, since when they are small 
enough they do not influence appreciably the dynamics, 
but mostly only the energies of the different modes in 
the steady state. In this case, the stimulated scattering 
constant R is irrelevant, as a change in its magnitude 
can be absorbed by the reservoir population ni^(x, t) and 
the condensate order parameter and results in 

the same equations with a renormalized pumping rate 
Pias' So, we can set R constant and vary only Pias- The 
threshold for condensation Pj^j^^ can be calculated numer- 
ically, since below P^j^^ the solution = is stable. 
The other parameters we used are m* = 6 x 10~^me, 

7c = C^ps)-\ = (?>mps)-K 

Our sample disorder potential is weak [12| , so the pro- 
flle of Fig. [5ja) is a resonable choice. The results depend 
on the choice of this potential, but there is no qualita- 
tive change as long as it has a local minimum and it is 
not perfectly symmetric. We also used a top-hat pump- 
ing proflle with radius of 8/im (Fig. [5jb)). Fig. [5fc-d) 
show the steady-state real space spectra calculated un- 
der two different pumping powers above threshold. Ex- 
cited states not conflned by the disorder potential spon- 
taneously appear inside the area pumped by the laser, in 
qualitative agreement with our experimental results. A 
detailed comparison with the experiment requires a more 
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